Decay of superfluid turbulence via Kelvin wave radiation and vortex reconnections. 
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The elementary processes involved in the decay of super- 
fluid turbulence in the limit of low temperature are studied 
by numerical simulations of vortex ring collisions. We find 
that small vortex rings produced by reconnections eventually 
annihilate in a collision where all their energy is converted 
into Kelvin waves and sound. We show that sound emission 
due to Kelvin waves also leads to a loss of vortex line length 
but this dissipation mechanism alone is too small to account 
for the experimentally observed decay. 

03.75.Fi, 67.40.Vs, 67.57.De 



One of the most challenging problems in fluid mechan- 
ics is the development of a complete hydrodynamical de- 
scription of turbulence. Progress is often made by study- 
ing fluids where the theoretical description can be simpli- 
fied. For example, low temperature superfluids such as 
liquid helium |l[] or the recently discovered dilute Bose- 
Einstein condensates Q are particularly attractive due to 
the dramatic reduction in viscosity. In a superfluid, the 
turbulent state consists of a network of interacting vor- 
tices known as a 'vortex tangle' In superfluid helium- 
4 the vortex tangle is strongly coupled to an interpene- 
trating normal fluid and the system shares many of the 
features of classical turbulence ||. However, in the limit 
of low temperature, the normal fluid is negligible, but 
experiments still indicate a temperature independent de- 
cay of the turbulence state [Q . In this case, the emission 
of sound either by vortex reconnections or vortex motion 
is the only active dissipation mechanism. Similarly, the 
time scale for the formation of vortex arrays in dilute 
Bose-Einstein condensates is also temperature indepen- 
dent suggesting that vortex-sound interactions may also 
be important in this relaxation process (^] . In an earlier 
work we showed that significant sound energy is released 
during vortex reconnections However, at the typical 
vortex line densities of superfluid hclium-4 experiments 
[4|, or in the regular structures found in a vortex lattice 
[5| one might expect that vortex reconnections are rel- 
atively infrequent, and that the continuous emission of 
sound by accelerating vortices may be the dominant dis- 
sipation mechanism. Unfortunately there is a scarcity of 
quantitative predictions of the significance of vortex mo- 
tion as a dissipation mechanism in superfluid turbulence. 
Analytical results for the sound radiated by moving vor- 
tices exist only for simple cases such as a co-rotating pair 
0. Furthermore, conventional numerical simulations of 
superfluid turbulence based on vortex filaments governed 



by incompressible Euler dynamics (the Biot-Savart law) 
H are unable to describe sound emission. 

An elegant model of quantum fluid mechanics is pro- 
vided by the Gross-Pitaevskii (GP) equation. The GP 
model represents an extension of Euler fluid mechanics to 
include the quantisation of circulation, vortex core struc- 
ture and vortex-sound interactions. Although the GP 
model does not accurately represent the physics of super- 
fluid helium-4, it has been shown to provide qualitative 
and sometimes quantitative insight into the critical veloc- 
ity for vortex nucleation fls^ JT^l, vortex line reconnections 
pT| , vortex ring collisions |12|, the decay of superfluid 
turbulence |jl~3f , and sound emission due to vortex recon- 
nections M. 

To study sound emission due to vortex motion we ex- 
cite Kelvin waves on a vortex ring and measure the length 
of the ring as a function of time. The Kelvin waves are 
produced by colliding two or more vortex rings. The ini- 
tial state is constructed from stationary vortex ring so- 
lutions of the uniform flow equation found by Newton's 
method Ejj. The desired configuration is obtained by 
multiplying the individual vortex ring states. This ini- 
tial state is then evolved according to the dimensionless 
GP equation, 



(1) 



using a semi-implicit Crank-Nicholson algorithm. In di- 
mensionless units, distance and velocity are measured in 
terms of the healing length, £, and the sound speed, c, 
respectively. In addition, the asymptotic number den- 
sity, no, is rescaled to unity. The computation box with 
volume, V — (50) 3 , is divided into 10 6 grid points with 
a spacing of 0.5. A grid spacing of 0.25 was also used to 
test the accuracy of the numerical methods. The time 
step is 0.02 and a typical simulation is run for 1.5 x 10 D 
steps. Simulations have been performed with two, three 
and four vortex rings with radii ranging from 2.86 to 
18.1 healing lengths. To convert the dimensionless units 
into values applicable to superfluid hclium-4, we take the 
number density as no = 2.18 x 10 28 m~ 3 , the quantum 
of circulation as n — h/m = 9.92 x 10~ 8 m 2 s _1 , and 
the healing length as £,/V2 = 0.128 nm |]|. This gives 
a time unit 27t£ 2 /k = 2 ps, and therefore for superfluid 
hclium-4 our simulations would correspond to a real time 
of 6 ns. Whereas for a sodium vapour condensate with 
£ ~ 0.2 /im [|J, the time unit is ~ 15 ns giving a simula- 
tion time of approximately 45 fis. 

In the first example we consider a collision between a 
large and a small vortex ring. Such collisions are impor- 



1 



tant in vortex tangle dynamics because small vortex rings 
arc often produced in collisions between larger structures 
|Q. A sequence of density isosurface plots illustrating 
the collision are shown in Fig. 1 Q Both vortex rings 
propagate in the positive x direction with their propaga- 
tion axes offset by the radius of the larger ring, i?i = 18.1. 
The collision occurs at t = 110, and leads to the de- 
struction of the small vortex ring. The vortex energy 
is converted into a sound pulse and Kelvin waves. The 
sound pulse propagates away from the large vortex ring 
and appears as a density minimum at t = 120 in Fig. 1. 
The Kelvin waves appear as two kinks which propagate 
around the vortex ring in opposite directions. 
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FIG. 1. Sequence of density isosurfaces (\ip\ 2 = 0.35) illus- 
trating a vortex ring collision with ring radii Ri = 18.1 and 
i?2 = 2.86. A top view is shown below each frame. The rings 
collide at t — 110. The reconnection produces a sound pulse 
which appears as a small dot in the t = 120 frame. 

In Fig. 2 we plot the vortex line length, £, as a function 
of time. The vortex line length is evaluated numerically 
by searching for phase singularities, estimating the length 
of vortex line within a grid cube and summing over all 
cubes. We have checked the accuracy of the method by 
performing the same calculation with a grid spacing of 



0.25 and 0.5. The two calculations agree to within a 
fraction of the healing length. During the reconnection 
at t = 110, the vortex line length is stretched from its 
initial value of 2w(Ri + R2) ~ 130, and then decreases 
dramatically due to the conversion of energy into sound 
||. After the collision the small vortex ring disappears 
completely while the mean radius of the large vortex ring 
is slightly larger, R[ = 18.8. 

The vortex line length oscillates due to the motion of 
the Kelvin waves around the vortex ring. By fitting the 
oscillations to a sine wave (shown bold in Fig. 2) we find 
that the period is 134 ± 1. This corresponds to a Kelvin 
wave velocity of 0.44 times the speed of sound. The dis- 
persion relationship for Kelvin waves is |ll| 
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where in dimensionless units, k — 2n and the core size 
a = 1/ V2, therefore the group velocity of a Kelvin wave 
with center wavelength A is 



2 In 



V2Y 



(3) 



A velocity of 0.44 corresponding to a Kelvin wavelength 
of 5.3 healing lengths, which agrees with a measure- 
ment made from Fig. 1. In this intermediate wave- 
length region, the Kelvin wave group velocity ([|) is only 
weakly dependent on wavelength (0.4 < v g < 0.63 for 
5 < A < 35). Consequently, the dispersion of the Kelvin 
wave packet is relatively slow. 

During the Kelvin wave oscillations the mean vortex 
line length gradually decreases. It is interesting to com- 
pare the observed length decrease with that expected for 
a Kelvin wave cascade, where energy is transfered to 
shorter and shorter wavelengths with a cut-off below a 
critical wavelength. In this case the vortex line density 
L = £/V, can be described by an equation of the form 
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where k = 2tt in dimensionless units, and \2 is a dimen- 
sionless coefficient. A linear fit to 1/L using the data 
from Fig. 2 for times between t = 500 and t = 2900 gives 
a coefficient, \2 ~ 0.006 ± 0.001. Both superfluid tur- 
bulence experiments |j"7| and numerical vortex dynamics 
simulations where an energy cut-off is assumed || suggest 
a coefficient, \2 ~ 0.3, in the limit of low temperature. 
This suggests that Kelvin wave radiation alone is not suf- 
ficient to account for the decay of superfluid turbulence. 



The o riginal animated movies for f igures 1 and 4 can be 
found at http://massey.dur.ac.uk/ml/. 
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FIG. 2. The vortex line length, £, as a function of time for 
the collision shown in Fig. 1. During the reconnection there 
is a sudden loss of line length due to the emission of a sound 
pulse. Subsequently the length decays due to Kelvin wave 
radiation. The bold line is a fit to a function of the form 
£ = £o[l + asin(2TYt/T + 0)]/(l + xrfot/V), where the period 
of the oscillations, T = 134±1, corresponds to the time for the 
Kelvin wave to go half way around the vortex ring. A linear 
fit to 1/L against t between t = 500 and t = 2900 determines 
the coefficient in Eq. (§), X 2 ~ 0.006 ±0.001. This data is for 
a grid spacing of 0.25 healing lengths. 

To test the generality of this result we repeated the 
calculation for different initial conditions. First, a colli- 
sion involving initial vortex rings with radii R\ = 18.1 
and i?2 = 5-74 produces the decay shown in Fig. 3. In 
this case the smaller vortex ring is not completely de- 
stroyed in the initial collision. Instead after making re- 
peated reconnections with the main vortex ring between 
t = 400 and t = 1250, it annihilates at t ~ 1250 leaving 
only the larger ring. The large vortex ring then decays 
by Kelvin wave radiation at a rate consistent with the 
results in Fig. 2. Second, collisions involving three vor- 
tex rings reproduce the same characteristic decay with 
large drops in the vortex line length during reconnection 
events interspersed with gently sloping regions of Kelvin 
wave radiation. All the simulations show that the Kelvin 
wave decay coefficient is independent of the initial condi- 
tion. This can be explained by the fact that the radiation 
spectrum is dominated by short wavelength Kelvin waves 




The length of vortex line destroyed during a recon- 
nection is strongly dependent on the reconnection angle 
. Consequently, to model a completely disordered vor- 
tex tangle we need to simulate a system which samples 
a complete range of reconnection angles. The simplest 
case where this is realised is four initial vortex rings, 
propagating inwards, as in the first frame of Fig. 4. 
The propagation axes are offset by one healing length 
to break the symmetry (the symmetric case completely 
annihilates in the first few hundred time units). A sim- 



ilar system, although with a very different length scale, 
has been studied using a classical vortex filament calcu- 
lation |l9| . The decay of the vortex line length is shown 
in Fig. 5. The initial vortex line density is much larger 
than in the previous examples, therefore there are a larger 
number of reconnections with a broad distribution of re- 
connection angles. A sequence of reconnections before 
t ~ 800 results in the formation of two large vortex rings 
which avoid each other until t ~ 2600. Again the decay 
coefficient for Kelvin wave radiation is consistent with 
Fig. 2. The collision at t ~ 2600 sets up another phase 
of repeated reconnections producing smaller vortex rings 
which annihilate in collisions with larger structures. This 
scenario of vortex tangle decay involving the production 
of small vortex rings is also observed in the classical vor- 
tex filament calculations of Tsubota et al. Q . Our main 
conclusion that Kelvin wave radiation alone cannot ac- 
count for the observed decay of superfluid turbulence is 
illustrated by plotting a decay characterised by \2 = 0.3 
in Fig. 5. 
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FIG. 3. The vortex line length as a function of time for 
a collision between vortex rings with radii Ri = 18.1 and 
i?2 = 5.74. The initial increase in the vortex line length is due 
to stretching during the first collision. Subsequently there are 
a large number of reconnections which result in the production 
of smaller rings. After t = 1100 a sequence of reconnections 
leaves only one vortex ring which then decays by Kelvin wave 
radiation. The decay between t = 1500 and t = 2500 is 
consistent with the coefficient found previously, \2 ~ 0.006. 

In summary, we have studied the basic processes in- 
volved in the decay of superfluid turbulence by numerical 
simulation of vortex ring collisions. At high vortex line 
densities our simulations suggest that the primary mech- 
anism for the loss of vortex line length is sound emission 
during reconnection events. The magnitude of the loss 
due to reconnections is in stark contrast to the relatively 
slow decay due to sound radiation from Kelvin waves. 
We conclude that Kelvin wave radiation alone is not suf- 
ficient to account for the decay of superfluid turbulence 
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but that Kelvin waves plus vortex reconnections could 
be. 
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FIG. 4. Sequence of density isosurfaces (IV"! 2 = 0.35) illus- 
trating a collision involving four converging vortex rings each 
with radius R — 18.1. A large number of reconnections result 
in the generation of small structures which annihilate into 
sound energy producing only two vortex rings by t — 740. 
Note that because of the periodic boundary conditions the 
loop structures connect on opposite faces of the numerical 
box. 



450 




1000 2000 
Time 



3000 



FIG. 5. The decay of vortex line length (thin line) for the 
collision involving four converging vortex rings shown in Fig. 
4. For t < 800 there are a large number of reconnections re- 
sulting in a dramatic decrease in length. After t ~ 800 only 
two rings remain and there are no more reconnections until 
t = 2400. In this region the loss of vortex line length fol- 
lows the prediction for Kelvin waves with a decay coefficient 
consistent with the data from Figs. 2 and 3. A decay char- 
acterised by X2 = 0.3 is plotted (thick line) to illustrate the 
difference between this case and the much slower decay rate 
due to Kelvin wave radiation only. 
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